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Abstract 

The extrapolation of nucleon magnetic form factors calculated within lattice QCD is investigated 
within a framework based upon heavy baryon chiral effective-field theory. All one-loop graphs are 
considered at arbitrary momentum transfer and all octet and decuplet baryons are included in 
the intermediate states. Finite range regularisation is applied to improve the convergence in the 
quark-mass expansion. At each value of the momentum transfer (Q 2 ), a separate extrapolation to 
the physical pion mass is carried out as a function of m, alone. Because of the large values of Q 2 
involved, the role of the pion form factor in the standard pion-loop integrals is also investigated. 
The resulting values of the form factors at the physical pion mass are compared with experimental 
data as a function of Q 2 and demonstrate the utility and accuracy of the chiral extrapolation 
methods presented herein. 

PACS numbers: 13.40.-f; 21.10.Ky; 12.39.Fe; ll.10.Gh 
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I. INTRODUCTION 



The study of the electromagnetic properties of the nucleon is of great importance in 
understanding the structure of baryons — see Refs. [l|, H, H, 0] for recent reviews. The most 
rigorous approach to low-energy phenomena in QCD is via numerical simulations in lattice 
gauge theory and many physical quantities, such as baryon masses, magnetic moments, etc. 
]!, @, 0, S H, E3, IH have been investigated within lattice QCD. Because of computing 
limitations, most of those quantities are simulated with large quark (tt) masses and an 
extrapolation of lattice results to the physical 7r mass is needed. Early lattice extrapolations 
considered simple polynomial functions of 7r mass. However, it is now widely acknowledged 
that the chiral non-analytic behavior predicted by chiral perturbation theory (x?T) must 



13|, ]1J, [15|, 116 



be incorporated in any quark mass extrapolation function [12|, 

%PT has been a very useful approach to the study of low momentum processes involving 
mesons and baryons and has been used in various studies of baryon structure. It is based 
on an effective Lagrangian constructed in a systematic way and consistent with all the 
symmetries of QCD. The first systematic discussion for the two flavor sector, i.e. the pion- 
nucleon system, of how to implement the ideas of chiral power counting 17], was performed 
in Ref. 18]. However, treating the nucleons as relativistic Dirac fields does not allow for 




a one-to-one correspondence between the expansion in small momenta and quark mass on 
the one hand and pion loops on the other. As pointed out in Ref. [HI], this shortcoming 
can be overcome if one makes use of methods borrowed from heavy quark effective field 
theory (HQEFT), namely to consider the baryons as extremely heavy, static sources. The 
relativistic or heavy baryon chiral perturbation theory has been applied to study a ran ge o f 
hadron properties in QCD, including nucleon magnetic moments and charge radii [2( 
the nucleon sigma commutator [22|, 23, 24| and moments of structure functions (l6l.]2 

Historically, most formulations of xPT are based on dimensional or infrared regulari- 
sation. However, the physical predictions of effective field theory must be regularisation 
scheme independent, such that other schemes are possible and may provide advantages 
over the traditional approach. Indeed, Donoghue et al. [29] have already reported the im- 
proved convergence of properties of effective theory formulated with what they called a 
"long-distance regulator". With the most detailed studies being on the extrapolation of 
the nucleon mass, it has been shown that the use of finite range regularisation (FRR) en- 
ables the most systematically accurate connection of ;\RT and lattice simulation results 
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23j, |26|, |27|, [ 

The FRR ;\RT has been applied to the extrapolation of proton magnetic moment with 



the leading non-analytic contribution of pions [3l|. It was found that the smooth behavior 



of the lattice data, together with the series truncations of the FRR expansion indicate that 
although higher order terms of DR can be individually large they effectively sum to zero in 
the region of interest. FRR x?T provides a resummation of the chiral expansion that ensures 
that the slow variation of magnetic moments observed in lattice QCD arises naturally in the 
FRR expansion. It was also predicted that the quenched and physical magnetic moments 
are in good agreement over a large range of pion mass, especially at large m n (3l| . 

In this paper, we will extrapolate the proton and neutron magnetic moments, as well as 
the form factors at finite momentum transfer, within a framework based upon heavy baryon 
chiral perturbation theory. Many methods have been used to compute the form factors at 
the physical value of the pion mass, including the early relativistic approach [l9j], heavy 
baryon chiral perturbation theory [32j], the so-called small scale expansion [33], relativistic 
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chiral perturbation theory [21|, |34|, etc.. The spectral functions of the form factors have 
also been investigated by calculating the imaginary parts of the form factors [35l . S. The 
disappointing observation was that a satisfactory description of the electromagnetic form 
factors was achieved only up to Q 2 ~ 0.1 GeV 2 . In order to improve this situation it 
is natural to consider higher order terms in chiral perturbation theory but eventually the 
series must diverge. In fact, effectively resumming the series by including vector meson 
degrees of freedom led to a satisfactory description of the electromagnetic form factors up 
to Q 2 
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It is of interest to calculate the form factors at relatively large 
momentum transfer and pion mass because there are many experimental data in this region 
and it is a priority for lattice QCD to understand that data. 

Because the values of the momentum transfer are quite large, it is not possible to make 
a systematic expansion in both Q 2 and m n . Instead, we extrapolate as a function of m n 
at each separate value of Q 2 . All the one loop contributions, including baryon octet and 
decuplet intermediate states, are considered. The quenched lattice data at large quark 
mass are used in the extrapolation, using the finding that the difference between quenched 
and full QCD data is usually quite small at large values of the pion mass. While this is a 
reasonable approach until high quality full QCD data is available (for first dynamical studies 



see Refs. [38l . l39l . |40j), it does mean that in comparing with experiment we must remember 
that there is an unknown systematic error associated with the use of quenched data. So 
too, we have not had lattice data available which would permit an extrapolation to the 
continuum (a — > 0) and infinite volume limits. In spite of all these caveats the results of 
this exploratory study are really very promising. 



II. CHIRAL PERTURBATION THEORY 



There are many papers which deal with heavy baryon chiral perturbation theory - for 
details see, for example, Refs. [4l|, 42,0, 44j. For completeness, we briefly introduce the 



formalism in this section. In the heavy baryon chiral perturbation theory, the lowest chiral 
Lagrangian for the baryon-meson interaction which will be used in the calculation of the 
nucleon magnetic moments, including the octet and decuplet baryons, is expressed as 

C v = iTrB v (v ■ V)B V + 2DTrB v S%{A^ B v } + 2FTrB v S^[A IM B v \ 

-if£(v ■ V)T vfl + C{f v ^A^B v + B v ApT£), (1) 

where is the covariant spin-operator defined as 

SS = ^7 5 ^X- (2) 

Here, v v is the nucleon four velocity (in the rest frame, we have v v = (1,0)). D, F and 
C are the coupling constants. The chiral covariant derivative is written as D^B V = 
d^B v + [V^,-Bt,]. The pseudoscalar meson octet couples to the baryon field through the 
vector and axial vector combinations 

v, = -(C^ + C^C), A, = i(C9 M c f - C%0, (3) 

where 

C = e**/ f , f = 93 MeV. (4) 
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The matrix of pseudoscalar fields is expressed as 

775*° + > ^ + 




7?° 



B v and T^ 4 are the velocity dependent new fields which are related to the original baryon 
octet and decuplet fields B and T M by 

B v (x) = e imN ^ x "B(x), (6) 

T*{x) = e imNVv ^T»(x). (7) 

In the chiral 577(3) limit, the octet baryons will have the same mass tub- In our calculation, 
we use the physical masses for baryon octets and decuplets. The explicit form of the baryon 
octet is written as 

B= E- -^° + 7E A n ■ (8) 

For the baryon decuplets, there are three indices, defined as 

T m = A + + T 112 = -J=A + , T 122 = i=A , (9) 
T 222 = A-, T 113 = ^S*- + , T 123 = -^=£*' , 

T223 — 7y^*' ' ^ 133 = "yf ' = '*' ' "^ 233 = ^/f " ' ' "^ 333 = ^ ' 

The octet, decuplet and octet-decuplet transition magnetic moment operators are needed 
in the one loop calculation of nucleon magnetic form factors. The baryon octet magnetic 
Lagrangian is written as: 

C = 4^ ^ dTtB ^ U Bv ) + ^B v a^ [F+ , B v ]) , (10) 

where 

Fi = \ (C^QC + CF^QC) ■ (ii) 

Q is the charge matrix Q =diag{2/3, —1/3, —1/3}. At the lowest order, the Lagrangian will 
generate the following nucleon magnetic moments: 

1 2 

Hp= ^HD + HF, Hn = —^HD- (12) 

The decuplet magnetic moment operator is expressed as 

g _ 

^ = -i—VcqijkT^uTvjuF^, (13) 
rriN 
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where q^ and qijkHc are the charge and magnetic moment of the decuplet baryon T^. The 
transition magnetic operator is 

C = ^/tF, v {e ljk Q]Bi m S^ klm + e^Qlf^S^B™) . (14) 



In Ref. [20|, the authors used fi u , Ha and fi s instead of the fic an d fir- F° r the particular 
choice, fi s = fid = — \fi u , one finds the following relationship: 

3 2 

A*d = = 3^0, fic = Hd, fMr = -4fi D - (15) 

In our numerical calculations, the above formulas are used and therefore all baryon magnetic 
moments are related to one parameter, /ip. 

In the heavy baryon formalism, the propagators of the octet or decuplet baryon, j, are 
expressed as 

v . k -5i N + is v ■ k - 5i N + is ' { ' 

where P^ u is v^v v — g^ u — (4/3)5^5^. 5 ab = mb — m a is the mass difference of between the 
two baryons. The propagator of meson j (j = n, K, rj) is the usual free propagator, i.e.: 

(17) 



k 2 - Mf + ie' 



III. NUCLEON MAGNETIC MOMENTS 

In the heavy baryon formalism, the nucleon form factors are defined as: 

< B{p')\J,\B{p) >= u(p') \v,G E {Q 2 ) + te ^ VaS " qU G M (Q 2 ))u(p), (18) 

where q = p' — p and Q 2 = —q 2 . According to the Lagrangian, the one loop Feynman 
diagrams which contribute to the nucleon magnetic moments are plotted in Fig. 1. The 
contributions to nucleon magnetic form factors of Fig. la are expressed as 

(la) _ m N (D + F) 2 NN m N (D + 3F) 2 /^ A + 3m N (D - F) 2 !^ 

~ 87T3/2 ^ + 487T3/2 > ^ 

m N (D + F) 2 

tNN m N (D - F) 2 



The integration I?f is expressed as 



8tT»£ - 87r3/ 2 

3 • A 



V — j dk — p , (21) 
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FIG. 1: The one loop Feynman diagrams for the nucleon magnetic moments. The solid, thick 
solid, dash and dotted lines are for the octet baryons, decuplet baryons, pseudoscalar mesons, and 
photons, respectively. 
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where 



Af = uj 3 (k + -q/2)uj(k - -q/2)(uj(k + "9 /2) + S af3 ) 

{O0j(k - -q/2) + 5 Q/3 )(c^(~fc + ~q/2) + Uj(k - ~q/2)). 



(22) 



ujj( k 



m 2 + k 2 is the energy of the meson j. In our calculation we use the finite range 

regularisation and u{ k ) is the ultra-violet regulator. This diagram is studied in the previous 
paper [3(| which gives the leading analytic term to the magnetic moments. The first terms 
in Eqs. (fl9j) and (120]) come from the ir meson cloud contribution. The second terms come 
from the K meson cloud contribution. Fig. lb is the same as Fig. la but the intermediate 
states are decuplet baryons. Their contributions to the magnetic form factors are expressed 
as 

m N C NA m N C fNTi * 



G 



P (ib) 



M 



G 



n(lfe) 
M 



367T3/2 

m N C 2 



NA 



m N C 2 



NT.* 



The contributions to the form factors from Fig. lc are expressed as 



(24) 



G 



,p(lc) 
M 



(D + F) 2 ([iEi — up) 

tNN 
~ 1 2tt 



1 



[D 



1927T3/2 



f)(f + 2f)Ml AS 



1927T3/2 



NT, 



D 

3 



(D- F) 2 (2^ + ^)1^ 

- F) 2 (hd + 3fi F ) NN 
192vr 3 /2 2^ , 



D 

3 



NA 



(25) 



G 



n(lc) 
M 



+ ( 



(D + F) 2 ^ F 
96vr 3 /2 

,2D 



T NN 
1 2tt 



1927T3/2 



(£> - F)*(jjl d - 2n F )I™ - (=- + F)Vd/^ 



2F)(D 



3 



where 



T a(3 
'2j 



d k 



967T3/2 



F ) 2 / i D jNN 



l 2 V 



10 A 



L 5j 



d k 



k){u j {k) + 
k 2 uCk) 2 



(26) 



(27) 



(28) 



The magnetic moments of baryons in the chiral limit, expressed in terms of fj,£> and fip, 
are used in the one loop calculations. However, we have taken the mass difference of the 
octet baryons into account. If the masses of the octet baryons are taken to be degenerate, 



then the coefficients in front of the integrals will be the same as in the paper of Ref. [41 
The contributions to the form factors of Fig. Id are expressed as 



G 



P (id) 

M 



hC 2 Hc T NA | 

i 27T "' 



G 



n(ld) 
M 



162tt 3 /2 

5C 2 /i C 
648vr 3 /2 



T NA 
1 2tt 



5C 2 yc 
1296tt3/2 j 2X 

5C 2 /i C 
1296vr 3 /2 



T NT* 



tNT* 
1 2K ■ 



(29) 



(30) 



7 



Fig. le and Fig. If give the following contributions to the form factors: 

r p(ie+if) _ (D + F)Cfx T NA 5(D- F)CjPr .jvss* , i£ + 3F)C/ir wvas* m, x 

r n(le+lf) _ (D + F)CfiT T NA , (-P ~ -F)C/X T ^j,* _ (-D + 3F)Cfl T NA %* 

Um ~ 1087r»/? 3 " 864tt3/2 5K 864tt3/2 , W 
where 

I 3 ? = f dt _> k2u{k J . (33) 

Fig. lg comes from the second order expansion of Lagrangian (fTUl) . The contributions to 
the magnetic form factors are expressed as 

327T3/2 hv-J^hK, (34) 



n n{lg) _ {^D + Vf) t {HO ~ Hf) t 

&M ~ 327T3/2 + 327r3/ 2 

where 



j 4j = / dk^H-. (36) 

The magnetic moment is defined as /i — Gm{Q 2 = 0). The total nucleon magnetic 
moments can be written as 

f, p (ml) =a p + aX + «K + £ ^(Q 2 = 0, mj), (37) 

AinW) = a n + a n 2 ml + a\m% + j^ G f k \Q 2 = 0,m^), (38) 
where (N = n, p) is expressed as 

J2GT\Q 2 = 0,ml = 0). (39) 



k=a 
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o N -r N 

k=a 



(?q (cq) is fj, p (fi n ) which is related to /xd and /if via Eq. (jl2"|h The residual series parameters, 
Oj, are determined by the best fit of the lattice data. 



IV. NUMERICAL RESULTS 

In the numerical calculations, the parameters are chosen as D = 0.76 and F = 0.50 
(gA = D + F = 1.26). The coupling constant C is chosen to be —1.2 which is the same 



as Ref. [41|. The renormalisation form factor u(k) can be monopole, dipole or Gaussian 



functions which give similar results [30j. In our calculations, the dipole function is used: 

u[k) = (ioW (40) 



s 



TABLE I: Residual series coefficients and nucleon form factors at zero momentum and 0.23 GeV 2 . 
The first two rows provide standard results whereas the final row includes the effect of the pion 
form factor. 



Q 2 


p 
% 


"0 


4 (GeV- 2 ) 


a n 2 (GeV- 2 ) 


a\ (GeV- 4 ) 


al (GeV" 4 ) 


M 







2.554 


-1.506 


-1.135 


0.420 


0.446 


-0.090 


2.73 ± 0.20 


-1.84 ± 0.19 


0.23 


1.617 


-0.932 


-0.411 


0.070 


0.144 


0.031 


1.70 ± 0.12 


-1.10 ± 0.11 


0.23 


1.652 


-0.968 


-0.499 


0.159 


0.201 


-0.027 


1.65 ± 0.10 


-1.06 ± 0.09 



with A = 0.8 GeV. The coefficients dp, aft, a>2, o>2, «4 and a n in Eqs. (1371) and (I38p are 
constrained by the quenched lattice data at large pion mass (m n > 500 MeV) where the 
quenched and physical values of the magnetic moments are expected to be close to each 
other @, [HI 



The K- and 77-meson masses have relationships with the pion mass according to 

2 ^ 2 , 2 I 1 2 1 I a 1 \ 

m K = + m K\ P hy - ^^Iphy, ( 41 ) 

m l = \ m l + m ?lphy - phy , (42) 

and enable a direct relationship between the meson dressings of the nucleon magnetic mo- 
ments and the pion mass. 

We begin by considering nucleon form factor results from the CSSM Lattice Collaboration 
To| . The proton magnetic moment fi p versus m 2 is shown in Fig. 2. Here, the last five 



lattice points at larger m 2 are used in the fit to avoid quenched chiral artifacts. The lines 
with label a and b+c+d+e+f+g correspond to the contributions of Fig. la and sum of the 
other diagrams, respectively. The residual series contribution, i.e. the contribution from 
a + c^m 2 + 04m 4 is also shown in the figure and labeled by "tree" . The near linear behavior 
of the residual series is a reflection of the excellent convergence of the residual expansion. 



The leading diagram (Fig. la) gives the dominant chiral behavior of the magnetic moment. 
At small pion mass, the proton magnetic moment decreases quickly with the increasing pion 
mass. At larger pion mass, the proton magnetic moment changes smoothly. At the physical 
point, = 0.139 GeV, the proton magnetic moment is 2.73 [In, close to the experimental 
value, 2.79 fiN- We emphasize again that the chiral curvature is dominated by Fig. la. 

The neutron magnetic moment, /i n , is studied in the same way. fi n versus m 2 is shown 
in Fig. 3. Again, only the five lattice points at larger pion mass are used in the fit to avoid 
quenched chiral artifacts. Similar to what was found in proton case, the leading diagram 
gives the dominant chiral curvature. The neutron magnetic moment increases quickly as 
ones moves from the chiral limit and becomes smooth at large pion mass. At the physical 
point, the neutron magnetic moment is —1.84 /1 at, compares favorably with the experimental 
value of — 1.91/iat. For both the proton and neutron, the rapid variation of the magnetic 
moments at = may reflect the fact that the nucleon magnetic radii diverge in the chiral 
limit. 
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FIG. 2: The proton magnetic moment [ICj] versus squared pion mass. The lines with label a and 
b+c+d+e+f+g correspond to the contributions of Fig. la and the sum of the other diagrams, 
respectively. 



In the above numerical calculations, we selected A to be 0.8 GeV. In Fig. 4 we show 
the nucleon magnetic moments versus A. The proton magnetic moment and the absolute 
value of neutron magnetic moment increase almost linearly with increasing A. In the range 
0.6 < A < 1.0 GeV, the proton (neutron) magnetic moment varies from 2.55 /ijy (—1-66 fi^) 
to 2.90 (—2.02 fi N ). When A is around 0.8 GeV, both the proton and neutron magnetic 
moments are in good agreement with the experimental values. 

Since the extrapolated values are A dependent (an indication that the fits lie outside the 
power-counting regime), the uncertainty of A will result in an additional source of error in 
the final result. Through a consideration of optimizing the convergence properties of the 
finite-range regularised expansion, we include the variation of A in the range 0.8 ± 0.2 GeV 
and add this uncertainty to the statistical uncertainties in quadrature. The extrapolated 
magnetic moments with corresponding error bars are listed in table I. 

In the chiral limit, the magnetic moments, c{J and Cq , are 3.41 and -2.53. These two values 
are close to the corresponding ones used in normal chiral perturbation theory. For example, 



in Ref. [341 ] . the corresponding values are 3.38 and -2.66. For the higher order terms, our low 
energy constants are much smaller, resulting in more convergent behavior. For example, our 
a? 2 and a 1 ^ are —1.14 and 0.42 which are much smaller than the corresponding values —6.80 
and 8.75 in Ref. [33 ] . 

We now proceed to extrapolate the nucleon magnetic form factors at finite Q 2 . At finite 
momentum transfer, we choose not to express the magnetic form factors in terms of Q 2 and 
Q as we did for the m n dependence of /iat. This is because the momentum dependence of 
form factors is close to the following assumption Gf /I (Q 2 ) = /ijv/(l + Q 2 /0.71GeV 2 ) 2 . The 
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FIG. 3: The neutron magnetic moment [13] versus squared pion mass. The lines with label a 
and b+c+d+e+f+g correspond to the contributions of Fig. la and the sum of the other diagrams, 
respectively. 
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FIG. 4: The proton and neutron magnetic moments versus the regulator scale A. 
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TABLE II: Residual series coefficients and nucleon form factors at various Q 2 . The first four lines 
report standard results while the latter four lines report results including the effect of a pion form 
factor. 



Q 2 (GeV 2 ) 


V 
a 


a" 


4 (GeV- 2 ) 


a n 2 (GeV- 2 ) 


rip 

M 




0.557 


1.042 


-0.638 


—0.024 


-0.00 


1.07 ± 0.17 


—0.71 ± 0.14 


1.08 


0.609 


-0.337 


0.015 


-0.04 


0.61 ± 0.13 


-0.37 ± 0.13 


1.14 


0.598 


-0.348 


0.052 


-0.04 


0.59 ± 0.11 


-0.37 ± 0.09 


2.28 


0.293 


-0.178 


0.035 


-0.03 


0.28 ± 0.09 


-0.18 ± 0.05 


0.557 


1.051 


-0.650 


-0.033 


0.01 


1.01 ± 0.15 


-0.66 ± 0.10 


1.08 


0.620 


-0.349 


0.008 


-0.03 


0.58 ± 0.12 


-0.34 ± 0.12 


1.14 


0.610 


-0.360 


0.044 


-0.04 


0.57 ± 0.10 


-0.35 ± 0.07 


2.28 


0.300 


-0.185 


0.032 


-0.02 


0.27 ± 0.09 


-0.17 ± 0.05 



high order terms in a Q 2 expansion are important and the truncation of momentum to some 
order, say fourth order, is not a good approximation. 

In our calculation, the same formulas as Eqs. (|3"7|) and (|38|) are used for the extrapolation 
of the magnetic form factors at each fixed finite value of Q 2 . The Q 2 dependence of nucleon 
magnetic form factors at tree level is included in the parameters do, a 2 and 04, as these 
parameters are constrained by the lattice results at finite Q 2 . 

In addition to the CSSM Lattice collaboration results l(j considered thus far, we also 
consider QCDSF lattice results at finite Q 2 [llj]. Large statistical uncertainties encountered 
at large Q 2 prevent one from constraining the m\ term and therefore we fit the QCDSF 
data using a residual series expansion up to and including order m 2 only. The coefficients 
together with the form factors at finite Q 2 are obtained by fitting the lattice results and 
are listed in Tables I and II. From the tables, one can see that a v Q decreases with increasing 
momentum, while increases with momentum, (N=p,n) are small indicating good 
convergence of the expansion. 

We plot the m n dependence of proton and neutron magnetic form factors at Q 2 = 0.23 
GeV 2 in Figs. 5 and 6 respectively. At small pion mass, the proton and neutron magnetic 
form factors do not change as quickly as in the case of zero momentum. However, the 
diagram of Fig. la still gives the dominant contribution to the curvature when the pion 
mass is small. From the figures, one can see that the expansion in powers of m n shows good 
convergence. At the physical pion mass, 6^(0.23 GeV 2 ) = (1.70 ± 0.12) /j, N and G^(0.23 
GeV 2 ) = —1.10 ± 0.11) fiN, which are both reasonable compared with experiment. 

In Figs. 7 and 8, we plot the proton and neutron magnetic form factors respectively. 
Results at Q 2 = 0.557, 1.08, 1.14 and 2.28 GeV 2 from the QCDSF collaboration Q are 
considered. From the figures, one can see that the lattice data do not vary smoothly as a 
function of the pion mass due to large statistical errors. As a consequence, the extrapolated 
magnetic form factors at the physical pion mass have relatively large error bars. Accurate 
lattice results are needed to better constrain the chiral expansion parameters and allow one 
to consider an m\ term in the residual expansion. 
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FIG. 5: The proton magnetic form factor 10( at Q 2 = 0.23 GeV 2 versus the squared pion mass. 
The lines with label a and b+c+d+e+f+g correspond to the contributions of Fig. la and the sum 
of the other diagrams respectively. 



To this point, we have not considered the possibility of an important role for the pion 
form factor in the calculation. We know that at large Q 2 , the pion form factor is much less 
than one and this will affect the meson cloud contribution to nucleon magnetic form factors. 
We list the pion electromagnetic form factor F n in Table III as provided in Ref. 45 . 



With these pion form factors, we repeat the chiral fit of the nucleon magnetic form factors. 
As an example, in Fig. 9, we show the result obtained for the proton magnetic form factor 
at Q 2 =0.23 GeV 2 . The dashed and solid lines are for the results with and without the pion 
form factor consideration, respectively. When the pion form factor is included, the leading 
diagram Fig. la provides less curvature. As a result, the total decreases from 1.70 ± 
0.12 fi N to 1.65 ± 0.10 /ijv- For the neutron, G^ increases from —1.10 ± 0.11 to —1.06 
± 0.09 /ijv. Though the pion form factor changes significantly at finite momentum, it does 
not affect the nucleon magnetic form factors very much. At large Q 2 , the pion form factor 
has a negligible effect on the nucleon form factors because the loop contribution itself is 
already very small. 

In Fig. 10, we show the extrapolated proton and neutron magnetic form factors versus Q 2 



TABLE III: Pion electromagnetic form factor, F w , at various Q 2 . 
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FIG. 6: The neutron magnetic form factor 10] at Q 2 = 0.23 GeV 2 versus squared pion mass. The 
lines with label a and b+c+d+e+f+g correspond to the contributions of Fig. la and the sum of 
the other diagrams respectively. 
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FIG. 7: The proton magnetic form factors [ll|] at Q 2 = 0.557, 1.08, 1.14 and 2.28 GeV 2 versus 



pion mass. 
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FIG. 8: The neutron magnetic form factors [11] at Q 2 = 0.557, 1.08, 1.14 and 2.28 GeV 2 
pion mass. 
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FIG. 9: The proton magnetic form factor 10|] at Q 2 = 0.23 GeV 2 versus pion mass. The dashed 
and solid lines are for the results with and without the pion form factor, respectively. 
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at the physical pion mass with the corresponding error bars. The hollow and solid square 
points are for the fits with and without the pion form factor, respectively. 

The solid lines in Fig. 10 are the empirical parametrization G^(Q 2 ) = Hn/(^ + 
Q 2 /0.71GeV 2 ) 2 . At both zero momentum and Q 2 = 0.23 GeV 2 , the extrapolated CSSM 
results are in good agreement with the experimental data. For the other values of Q 2 from 
the QCDSF collaboration, the extrapolated nucleon magnetic form factors are in reasonable 
agreement with the empirical parameterization. For the neutron the agreement is quite rea- 
sonable, while for the proton, although the extrapolation is consistently within one standard 
deviation of the empirical curve, the extrapolated values do appear to be systematically a 
little high. 

We should mention that we use just two parameters, ao and 02, to fit the lattice data (at 
each value of Q 2 ). We note that it would be very helpful to have more accurate lattice data 
over a range of lattice spacings and volumes in order to extrapolate to the infinite volume 
continuum limit and to be able to incorporate an a& term. One would also prefer to work 
with full QCD data rather than quenched data. Until these conditions are satisfied it is a 
little early to draw strong conclusions about the validity of the extrapolation in pion mass 
from a comparison with experimental data for the form factors. Indeed, that the current 
results lie within one standard deviation of the data at all values of Q 2 is really a very 
positive result at the present stage. We do emphasis that our results are based on the lowest 
order Lagrangian at one loop level in heavy baryon approximation. That must breakdown at 
high momentum transfer, however, that is precisely where, in the FRR treatment, the loops 
become naturally small - for the clear physical reason that high pion momenta (and high 
pion mass) are suppressed by the finite size of the source. The extrapolation error arising 
from the uncertainty in A included above is small when the momentum is high. Again, this 
is because at high momentum transfer, the loop contribution to the total magnetic form 
factor is small. 

V. SUMMARY 

We extrapolated state of the art lattice results for nucleon magnetic form factors in an 
extension of heavy baryon chiral perturbation theory. All one-loop graphs are considered 
at arbitrary momentum transfer and all octet and decuplet baryons are included in the 
intermediate states. 

Finite-range regularisation is used in the one loop calculation to improve the conver- 
gence of the chiral expansion. The residual series coefficients a^a^), °2 (°2) an d a \ (04) 
are obtained by fitting the lattice results at m T > 0.5 GeV, where quenched artifacts are 
anticipated to be small. 

The leading non-analytic diagram provides the dominant curvature for the m n depen- 
dence of magnetic moments. The sum of higher-order one-loop terms provide only a small 
correction to this curvature. The one-loop contributions show that the proton (neutron) 
magnetic moment decreases (increases) quickly with increasing pion mass in the small m n 
region. At larger pion masses, their contributions change slowly and smoothly. 

The magnetic form factors are also studied at large Q 2 where chiral nonanalytic behavior 
is suppressed. Here, the importance of the pion form factor is also examined. For Q 2 = 0.23 
GeV 2 , G P M = 1.70 ± 0.12 fi N and G n M = -1.10 ± 0.11 fi N . Upon including the pion form 
factor, these values will change to 1.65 ± 0.10 [i N and —1.06 ± 0.09 fi N indicating the effect 
is subtle. Although the pion form factor decreases quickly with increasing momentum, its 
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FIG. 10: The proton and neutron magnetic form factors at physical pion mass versus the momen- 
tum transfer, Q 2 . The hollow and solid square points are for the fits with and without the inclusion 
of the form factor of the pion, respectively. 



effect on nucleon form factors is not significant, as the loop integrals themselves are already 
small. 

The chirally extrapolated results of Fig. 10 compare favorably with experiment and 
demonstrate the utility and accuracy of the chiral extrapolation methods presented herein. 
With the more accurate lattice data and treatment, it is of interest to see whether the 
mismatch at large Q 2 will disappear or not. 



Acknowledgements 

P.W. thanks the Theory Group at Jefferson Lab for their kind hospitality. This work 
was supported by the Australian Research Council and by DOE contract DOE-AC05- 
06OR23177, under which Jefferson Science Associates operates Jefferson Lab. 



[1] H. Y. Gao, Int. J. Mod. Phys. E 12, 1 (2003) [Erratum-ibid. E 12, 567 (2003)] 
[arXiv:nucl-ex/0301002~] . 

[2] C. E. Hyde-Wright and K. de Jager, Ann. Rev. Nucl. Part. Sci. 54, 217 (2004) 
|arXiv:nucl-ex/050700T| . 



[3] J. Arrington, C. D. Roberts and J. M. Zanotti, arXiv:nucl-th/0611050 



[4] C. F. Perdrisat, V. Punjabi and M. Vanderhaeghen, arXiv:hep-ph/0612014 



17 



[5] D. B. Leinweber, W. Melnitchouk, D. G. Richards, A. G. Williams and J. M. Zanotti, Lect. 

Notes Phys. 663, 71 (2005) |[ arXiv:nucl-th/0406032] . 
[6] D. B. Leinweber, R. M. Woloshyn and T. Draper, Phys. Rev. D 43, 1659 (1991). 
[7] D. B. Leinweber et ai, Phys. Rev. Lett. 94, 212001 (2005) |arXiv:hep-lat /0406002 . 
[8] D. B. Leinweber et ai, Phys. Rev. Lett. 97, 022001 (2006) [arXiv:hep-lat /0601025 . 
[9] R. G. Edwards et al [LHPC Collaboration], Phys. Rev. Lett. 96, 052001 (2006) 

[arXiv:hep-lat /0510062] . 

[10] S. Boinepalli, D. B. Leinweber, A. G. Williams, J. M. Zanotti and J. B. Zhang, Phys. Rev. D 
74, 093005 (2006). 

[11] M. Gockeler et al [QCDSF Collaboration], Phys. Rev. D 71, 034508 (2005) 

[arXiv:hep-lat/03030l9~] . 
[12] D.B. Leinweber and T.D. Cohen, Phys. Rev D 47, 2147 (1993). 
[13] D.B. Leinweber, AW. Thomas and R.D. Young, Phys. Rev. Lett. 86, 5011 (2001). 
[14] T.R. Hemmert and W. Weise, Eur. Phys. J. A 15, 487 (2002). 

[15] D. B. Leinweber, A. W. Thomas, K. Tsushima and S. V. Wright, Phys. Rev. D 61, 074502 

(2000) [arXiv:hep-lat /9906027| . 
[16] W. Detmold et al, Phys. Rev. Lett. 87, 172001 (2001) [arXiv:hep-lat /0103006| . 
[17] S. Weinberg, Physica 96A (1979) 327. 

[18] J. Gasser, M. E. Sainio and A. Svarc, Nucl. Phys. B 307, 779 (1988). 
[19] E. Jenkins and A. V. Manohar, Phys. Lett. B 255, 558 (1991). 

[20] P. Ha and L. Durand, Phys. Rev. D bf 58, 093008 (1998); Phys. Rev. D 67, 073017 (2003). 
[21] B. Kubis and U. G. Meissner, Nucl. Phys. A 679, 698 (2001). 

[22] B. Borasoy and U. G. Meissner, Annals Phys. 254, 192 (1997) [arXiv:hep-ph/9607432| . 
[23] D. B. Leinweber, A. W. Thomas and R. D. Young, Phys. Rev. Lett. 92, 242002 (2004). 
[24] M. Procura, T. R. Hemmert and W. Weise, Phys. Rev. D 69, 034505 (2004). 
[25] T. R. Hemmert, M. Procura and W. Weise, Phys. Rev. D 68, 075009 (2003). 
[26] D.B. Leinweber, A. W. Thomas, K. Tsushima and S. V. Wright, Phys. Rev. D 61, 074502 
(2000). 

[27] W. Armour, C. R. Allton, D. B. Leinweber, A. W. Thomas and R. D. Young, J. Phys. G 32, 

971 (2006) | arXiv:hep-lat/0510078| . 
[28] C. R. Allton, W. Armour, D. B. Leinweber, A. W. Thomas and R. D. Young, Phys. Lett. B 

628 (2005) 125 | arXiv:hep-lat /0504022| . 
[29] J. F. Donoghue, B. R. Holstein and B. Borasoy, Phys. Rev. D 59, 036002 (1999). 
[30] R. D. Young, D. B. Leinweber and A. W. Thomas, Prog. Nucl. Phys. 50, 399 (2003). 
[31] R. D. Young, D. B. Leinweber and A. W. Thomas, Phys. Rev. D 71, 014001 (2005) 

|arXiv:hep-lat/0406001~] . 
[32] V. Bernard, N. Kaiser, J. Kambor and U. G. Meissner, Nucl. Phys. B 388, 315 (1992). 
[33] V. Bernard, H. W. Fearing, T. R. Hemmert and U. G. Meissner, Nucl. Phys. A 635, 121 

(1998) [Erratum-ibid. A 642, 563 (1998)] [arXiv:hep-ph/9801297| . 
[34] T. Fuchs, J. Gegelia and S. Scherer, J. Phys. G 30, 1407 (2004) |arXiv:nucl-th/0305070| . 
[35] V. Bernard, N. Kaiser and U. G. Meissner, Nucl. Phys. A 611, 429 (1996) 
[36] N. Kaiser, Phys. Rev. C 68, 025202 (2003) [arXiv:nucl-th/0302072] . 

[37] M. R. Schindler, J. Gegelia and S. Scherer, Eur. Phys. J. A 26, 1 (2005) 
[arXiv:nucl-th/0509005| . 

[38] C. Alexandrou, G. Koutsou, J. W. Negele and A. Tsapalis, Phys. Rev. D 74, 034508 (2006) 
|arXiv:hep-lat/06050l7] . 



18 



[39] M. Gockeler et ai, |arXiv:he p-lat/060900T| 
[40] R. G. Edwards et ai, |arXiv:hep-lat/0610007| 

[41] E. Jenkins, M. Luke, A. V. Manohar and M. J. Savage, Phys. Lett. B 302, 482 (1993); 

Erratum-ibid. B 388, 866 (1996). 
[42] J. N. Labrenz and S. Sharpe, Phys. Rev. D 54, 4595 (1996). 
[43] L. Durand, P. Ha, and G. Jaczko, Phys. Rev. D 64, 014008 (2001). 
[44] B. C. Tiburzi, Phys. Rev. D 71, 054504 (2005) | arXiv:hep-lat /0412025] . 



[45] J. P. B. C. de Melo, J. S. Veiga, T. Frederico, E. Pace and G. Salme, hep-ph/0609212 



19 



